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ABSTRACT 

We study the propagation effects of radio waves in a pulsar magnetosphere, com- 
posed of relativistic electron-positron pair plasmas streaming along the magnetic field 
lines and corotating with the pulsar. We critically examine the various physical ef- 
fects that can potentially influence the observed wave intensity and polarization, in- 
cluding resonant cyclotron absorption, wave mode coupling due to pulsar rotation, 
wave propagation through quasi-tangential regions (where the photon ray is nearly 
parallel to the magnetic field) and mode circularization due to the difference in the 
electron/positron density/velocity distributions. We numerically integrate the transfer 
equations for wave polarization in the rotating magnetosphere, taking account of all 
the propagation effects in a self-consistent manner. For typical magnetospheric plasma 
parameters produced by pair cascade, we find that the observed radio intensity and 
polarization profiles can be strongly modified by the propagation effects. For relatively 
large impact parameter (the minimum angle between the magnetic dipole axis and the 
line of sight), the polarization angle profile is similar to the prediction from the Rotat- 
ing Vector Model, except for a phase shift and an appreciable circular polarization. For 
smaller impact parameter, the linear polarization position angle may exhibit a sudden 
90° jump due to the quasi-tangential propagation effect, accompanied by complex cir- 
cular polarization profile. Some applications of our results are discussed, including the 
origin of non-gaussion pulse profiles, the relationship between the position angle pro- 
file and circular polarization in conal-double pulsars, and the orthogonal polarization 
modes. 
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1 INTRODUCTION 

Pulsar radio emission is likely generated within a few hundred kilometers from the neutron star (NS) surface (e.g. Cordes 1978; 
Blaskiewicz et al. 1991; Kramer et al. 1997; Kijak & Gil 2003). A pulsar is surrounded by a magnetosphere filled with relativistic 
electron-positron pair plasmas (plus possibly a small amount of ions) within the light cylinder. When radio waves propagate 
though the magnetosphere, the total flux, polarization state and spectrum of the emission may be modified by propagation 
effects. Understanding the property of wave propagation in pulsar magnetospheres is necessary for the interpretation of various 
observations of pulsars. 

Radio emission from pulsars shows strong linear polarization. For some pulse components or even the whole pulse profiles 
it can be 100% percent polarized (e.g. Lyne & Manchester 1988; Gould & Lyne 1998; Weisberg et al. 1999, 2004; Han et al. 
2009). Linear polarization (LP) is closely related to magnetic field lines where the emission was generated. Based on the linear 
polarization position angle (PA) curve of Vela pulsar, the Rotating- Vector-Model (RVM) was suggested by Radhakrishnan 
& Cooke (1969). For some pulsars, especially the so-called conal-double type pulsars, RVM works very well (e.g., Mitra & 
Li 2004). However, the PA curves of most pulsars are much more complex and do not follow the simple RVM model. The 
deviation from the RVM model could be caused by the intrinsic emission mechanism (e.g., Blaskiewicz et al. 1991), which is 
highly uncertain (e.g., Lyubarsky 2008), and/or the propagation effect through the pulsar magnetosphere (see below). Also, 
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the PA curves or polarization observations of individual pulses show the orthogonal polarization modes (OPM) phenomenon, 
in which the polarization position angle exhibits a sudden ~ 90° jumps (e.g. Manchester et al. 1975; Backer et al. 1976; Cordes 
et al. 1978; Stinebring et al. 1984a, 1984b; Xilouris et al. 1995). It is not clear whether the OPM arises from the emission 
process (e.g. Luo & Melrose 2004) or the propagation effect (e.g. McKinnon & Stinebring 2000) . 

Another important observational feature of pulsar radio emission is the circular polarization (CP, e.g. Rankin 1983; 
Radhakrishnan & Rankin 1990; Han et al. 1998). Significant CPs have been observed in individual pulses of pulsars with 
mean values typically 20%-30%. Very high degrees of CP are occasionally observed from some components of pulsar profiles 
(e.g. Cognard et al. 1996; Han et al. 2009). Radhakrishnan & Rankin (1990) identified two main types of CP signature: 
antisymmetric type with sign reverse in the mid-pulse and symmetric type without sign change over whole profile. They 
concluded that the CP of the antisymmetric type is associated with the core emission and strongly correlated with the sense 
of rotation of the linear position angle. Han et al. (1998) showed that this correlation is not kept for a larger sample, and they 
found that for conal-double pulsars the sense of CP is correlated with the sense of PA curves. 

The diverse behaviours of pulsar polarization (including LP and CP) may require more than one mechanisms for proper 
explanations. First, they may be caused by an intrinsic mechanism in the emission region and/or process. For example, 
Randhakrishnan & Rankin (1990) suggested that geometrical effect to the pulsar beam from curvature radiation can naturally 
generate antisymmetric circular polarization for the core components. Gangadhara (1997) suggested that the observed circular 
polarization could be caused by the coherent superposition of two orthogonal modes emitted by positrons and electrons. Xu 
et al. (2000) interpreted the circular polarization by the superposition of coherent inverse Compton scattering. Kazbegi 
et al. (1991) suggested that cyclotron instability may be responsible for the circular polarization. Also, Luo & Melrose 
(2001) suggested that circular polarization can develop by cyclotron absorption when the distributions (especially the number 
densities) of the magnetospheric electrons and positrons are different. 

However, many observed characteristics of the pulsar radio emision are most likely dictated by the wave propagation in 
the magnetospheric plasma (see, e.g., Melrose 2003 and Lyubarsky 2008 for a review). A number of theoretical works have 
been devoted to study how magnetosphere propagation influences pulsar polarization observations. Whatever the emission 
mechanism, radio wave propagates in the plasma in the form of two orthogonally polarized normal modes. The polarization 
state of the wave evolves along the ray, following the direction of the local magnetic field, a process termed "adiabatic walking" 
(Cheng & Ruderman 1979). Cheng & Ruderman (1979) introduced two propagation effects: the wave mode coupling effect 
for pure pair plasma and the circularization effect (natural modes become circular polarized), both of which can generate 
circular polarization. Melrose (1979) and Allen & Melrose (1982) suggested that the separation of natural waves (because of 
different refractive indices) can cause the OPM phenomenon. Arons &: Barnard (1986) studied the wave dispersion relation 
and natural modes in the relativistic pair plasma. Lyubaskii & Petrova (1999) considered the natural modes in relativistic 
plasma with co-rotating velocity in the infinite magnetic field limit, and Petrova & Lyubarskii (2000) studied refraction and 
polarzation transfer in such a plasma. Luo & Melrose (2001) and Fussell et al. (2003) studied the cyclotron absorption of radio 
emission within pulsar magnetospheres. Petrova (2006) further studied the polarization transfer in pulsar magnetosphere and 
considered the wave mode coupling and cyclotron absorption effect. Johnston et al. (2005) suggested that the variation of 
circular polarization of PSR B1259— 63 during the elipse with its main-sequence companion is related to the wave propagation 
effect in the magnetosphere of the companion star. However, none of the previous studies have calculated the final polarization 
profiles with all of these propagation effects included in a self-consistent way within a single theoretical framework. It is often 
unclear which of the effects are most important, and if so, under what conditions. In this paper we attempt to combine all 
the propagation effects, evaluate their relative importance, and use numerical integration along the photon ray to study the 
influence of propagation effects on the final polarization states. 

This paper is organized as follows. In section 2, we present the geometrical model for our calculation and the general 
wave evolution equation in a magnetized plasma. In section 3, we give the expression of the dielectric tensor of a relativistic 
pair plasma characterizing the magnetosphere of a pulsar, and discuss the natural wave modes and their evolution. In section 
4, we study several important propagation effects separately: cyclotron absorption, wave mode coupling, circularization and 
the quasi-tangential propagation (see Wang & Lai 2009). In section 5, we present numerical calculations of the single photon 
evolution and the phase profiles of pulsar emission beam. Our results and possible applications are presented in section 6. 



2 GEOMETRY AND GENERAL WAVE EVOLUTION EQUATION 
2.1 Geometrical Model 

Consider a photon (radio wave) emitted at the initial position r; at time U (corresponding to the pulsar rotation phase ^i). 
Suppose the photon trajectory is a straight line along k (the wave vector). In a fixed XYZ frame with Z — k along the line 
of sight and f2 (the pulsar spin vector) in the XZ-p\&ne (Z x tl = svciC^Y , here C, is the angle between k and ft; see Fig. 1), 
the photon position after emission and the corresponding pulsar rotation phase are 
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Figure 1. Three frames used in this paper: 1) The fixed frame XYZ with Z |j fc, f2 in the XZ-plane and fc X fi = fishily. The 
direction of B in this frame is (8 b, 4>b)', 2) The instantaneous incrtial frame xyz with z = Z, B in the rrz-plane and fc X B = — sin 9 b y\ 
3) The instantaneous inertial frame x'y' z' with z' = B, k in the rr'z'-plane and fc X B = — sin 8bv' ■ 



r = n + sZ, (2.1) 

* = *i + Q{t - U) = *i + s/ric, (2.2) 

where s = c(t — U) is the distance from the emission point along the ray, and r\ c — c/f2 the radius of the light cylinder. The 
rotating magnetic field is given by 

B(s) = -V( M ■ r/r 3 ) = -ii + • r), (2.3) 

with 

(j,(s) = /i [(sin ("cos a — cos ("sin a cos $!)X — sin a sin ^Y + (cos ( cos a + sin ("sin a cos "t)^] , (2.4) 

where a is the inclination angle between fi and fi (see Fig. 1). Note that the impact angle x, which is the smallest angle 
between fc and /x, is given by \ — C ~ a - Thus, the polar angles of /x in XYZ frame, (# M , <^) M ), are given by 

— sin a sin VP 
sin £ cos a — cos £ sin a cos * ' 

The magnetic field at a given point along the ray is inclined at an angle 9b with respect to the line of sight, and make 
an azimuthal angle 4>b in the XY-plane such that: 

cos6> s (s) = tan<^) S (s) = (2.6) 

a tlx 



cos flj, = cos ( cos a + sin ( sin a cos W, tan</> (J = — — — — . (2.5) 



2.2 Wave Evolution Equations 

The wave equation for photon propagation takes the form 

2 

Vx (fi' 1 -V x E) = ^je- E, (2.7) 

where E is the electric field, and e, /x -1 are the dielectric and inverse permeability tensors, respectively. The inverse perme- 
ability is very close to unity when B <C Bq = 4.414 x 10 13 G (the critical QED field strength), and we set /x to be unity in the 
remainder of the paper. In practice, it is most convenient to calculate dielectric tensor in the x'y' ' z frame (where the z'-axis 
is along B, and k in the 2;'z'-plane, see Fig. 1). Once [e] x , y , z , (the matrix representation of the dielectric tensor in the x'y 1 V 
frame) is known, we can easily obtain [e] XYZ in the fixed XYZ frame through a coordinate transformation 

[e] XYZ = M[e] xlylz ,M T , (2.8) 

where the transformation matrix M is 
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-cos 0s sin 0s 
M — I —sin 0s —cos 0s 





- cos 6<s cos 0s sin 0s sin #s cos 0s 
= | — cos 6 b sin 0s —cos 0s sin #s sin 0s 
sin#s cos 6 b 



(2.9) 



and M T is the transpose matrix of M. 

Knowing e along the trajectory, we can use eq. (2.7) to derive the wave amplitude evolution equation. Let E = e lk ° s A, 
where k = ui/c. Assuming that |cL4/ds| -C fco|A| (geometric optics approximation), we obtain 
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(2.10) 



(2.11) 

The wave evolution equation (2.10) can be used to study the evolution of EM wave amplitude across the pulsar magnetosphere. 

We can also follow the evolution of the four Stokes parameters instead of the evolution of the wave amplitudes. The four 
Stokes parameters are defined by (in the fixed XY Z frame) 

I = AxAx+AyAy, 
Q = A X A X — A Y A Y , 
U = AxAy+AyAx, 

V = -i(AxAy-AyAx). (2.12) 
Combining with eq. (2.10), we obtain the evolution equations for the Stokes parameters: 
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(2.13) 
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Here the subscript "i" and "r" correspond to the real and imaginary part of each element. 

If we know the dielectric tensor along the ray, we can integrate eq. (2.10) from the emission point in the inner magneto- 
sphere to large distance where the plasma no longer affect the radiation (both intensity and polarization). We will calculate 
the dielectric tensor of a relativistic streaming pair plasma in the next section. 



3 WAVE MODES AND PROPAGATION IN A STREAMING PLASMA 

The magnetospheres of pulsars consist of relativistic electron-positron pair plasma streaming along magnetic field lines. The 
Lorentz factor 7 of the streaming motion and the plasma density iV are uncertain. For the open field line region of radio 
pulsars, pair cascade simulations generally give 7 ~ 10 2 — 10 4 and r\ = N/Nqj ~ 10 2 — 10 5 (e.g., Daugherty & Harding 
1982; Hibschman & Arons 2001; Medin & Lai 2009), while recent theoretical works suggest that the corona of magnetars 
consist of pair plasma with 7 up to 10 3 and 77 ~ 2 x 10 3 (R*/r) (where R, is the stellar radius; Thompson et al. 2002; 
Beloborodov & Thompson 2007). Here Ngj = (£IB) / \2ixec) is the Goldreich- Julian density. In this paper, we choose plasma 
density r\ = N/Ngj to be in the range of 100 - 1000, and the Lorentz factor 7 of the streaming motion to be 100 - 1000. 
We also consider a small asymmetry between positrons and electrons, i.e. AN/N 7^ and A7/7 7^ 0, where AJV, A7 are the 
differences in the number densities and Lorentz factors between electrons and positrons. 



3.1 Dielectric tensor 



The dielectric tensor e in the x'y'z' frame (with z' = B, k in the a/z'-plane and k x B = — sinSsy'; see Fig. 1) can be 
written as [see eqs. (2.11) - (2.13) and (2.19) of Wang & Lai 2007]: 
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£xj II 1 
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ty'x' 
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(3.14) 



ty'z' = < 
t-z'z' 

with 

/s, 11 = 

/s,12 = 

/s,ij = 
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^2 J fsMfs(ls)d'ys, 
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S 
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(3.15) 



«s7s 1 (1 + «7rad) 



(1 + I7rad) 2 - Us'fs 2 (1 - P 3 COS B ) 2 ' 

sign (q s ) ul /2 v s ~/~ 2 (1 - f3 a cos fla) -1 
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Here the subscript "s" specifies different species ("e" is for electron and "p" for positron), and /3 S , 7 S and / s (7s) are the 
velocity (divided by c), Lorentz factor and its distribution function. The dimenssionless parameters u 3 , v s are 

(3.17) 
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(3.18) 

Here N s is the number density of particles, N = N p + N e , lo c and u p \ are the cyclotron and plasma frequencies, which are 
given by 

^ = ujc_ = J_eB_ = 2 7Q5 x 1q9 QHz 

2w Z7T m e c 



fpl = 



47riVe 2 



= 8.960 x 10 3 iV 1/2 Hz = 2.370 ri L ' z B L 1 ^P' 1 



l/2 R l/2p-l/2 



GHz, 



(3.19) 
(3.20) 



2-7T 27T V 

where the magnetic field B12 = -B/(10 12 G), the pulsar spin period P lB — P/(ls), and the dimensionless density r\ — N/Nqj 
is measured in units of the Goldreich- Julian density, JVgj = Q.B/(2irec) ~ 7.0 x 10 10 £>i2-Pis cm -3 . The refractive index, 
n — ck/u, is generally very close to unity, so we always set n ~ 1 here. The radiative damping 



7rad 



4e 2 ^ c 
3mc 3 



(3.21) 



is important only near the cyclotron resonance [where 7s^(l — n/^cos^s) ~ uj c ] and can be neglected at other places. The 
function sign(q s ) equals to —1 for electrons and 1 for positrons. 

In this paper we focus on cold streaming plasmas, which means that both electrons and positrons in the streaming plasma 
have single 7 S or / s (7s) = <5(7s _ 7s,o)- Thus, we need not to integrate across 7 S when calculating each element of the dielectric 
tensor in eq. (3.15). 

When we consider the region r <C r cyc (the cyclotron resonance radius), we can take the infinite magnetic field limit, and 
the damping term can be neglected. In this case the dielectric tensor becomes very simple (e.g., Arons & Banard 1986) 



\^\x' y' z' 



1 
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1 + fr, 



(3.22) 
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with /„ = -wy- s (l - f3cos9 B )-- 

3.2 Wave evolution equation for single-7 plasma 

In this subsection we consider the polarization evolution equation for single 7 plasma, i.e. all electrons (positrons) have the 
same 7e(7p)- We assume that there is a small asymmetry between electrons and positrons in N 3 or 7s : AN/N <JC 1, where 
N — N p + N e , AN = N p — N e (usually AN/N is the reciprocal of the multiplicity of the cascade), and/or A7/7 <JC 1 [where 
7 = (-y p + 7e)/2, A7 = jp — 7e]. In this case the final matrix elements in the wave evolution equation (2.10) are 

oxx = -Fii(l + fe cos 2 4> B ) + F v sin 2 9 B cos 2 (f> B — + F v sin 2 9 B cos 2 <f> B , 

<?xy = Fufe + F v sin 2 B sin^s cos(j> B — iF\2 ~ F v sin 2 9 B sin^s cos<j> B — iF\2, 

pyx = Fn fg + F v sin 2 6 B sin 4> B cos (/) B + iFi2 ~ F v sin 2 6 B sin 4> B cos </) B + iFi2, 

oyy = F\\{\ + fe sin 2 </> B ) + F v sin 2 9 B sin 2 <j) B — Fu + F v sin 2 6 B sin 2 </> B , (3.23) 



with 



Ef 7 
/mi — ~ 1 - JTi a 'H~TT = Fll > r + ^n.*! 



F12 



EST^ sign(q s )v s y s ul /2 -/ s (l-f3 s cosb B ) 

^ l + 2j 7r ad-^7s" 2 l~/3 S COSe s r 2 



+ 2i7 rad -u s 7 s (1 -- /3 S cos6> B )~ 



F " = XI /s '" - -«7 _3 (l-/3cos^ s )- 



2 



/, = (cosg B ^Csing B ) 2 -l^- (1 ^ 2)2 . (3.24) 

In the deriving of eq. (3.23), we have assumed #57 3> 1 (which is valid for most places), so that fg ~ 0. 
Using eqs. (2.13) and (3.23), we can write the evolution equation of the four stokes parameters as 

^ = -fcoFu.i/ - k F 12 ,iV, 
as 

^ = -fc o F 11 , l Q + fc o F 12 , r [/+yF I) sin 2 6l s sin20 s \/, 

= -fc Fl2.r<3 - fcoFn.it/ - ^F,, Sin 2 6»s COS 20BV, 

as 2 

^ = -fcoF^,*/- yF„sin 2 fe(Qsin2<^B - U cos 2(f) b) - k Fu,iV. (3.25) 

Here fco = c/u, the subscript "r" and "i" specify the real and imaginary parts. Equation (3.25) is useful for understanding 
the different kinds of propagation effects on the polarization evolution (see section 4) . 

3.3 Wave modes 

Using the electric displacement D = e-E in the Maxwell equations, we obtain the equation for plane waves with E oc e '0* , -" t ) 

[etj + n 2 (kikj - Sij)]Ej = 0, (3.26) 

where n = ck/to is the refractive index and k — k/k. In the coordinate system xyz with k along the 2-axis and B in the 
xz-plane (see Fig. 1), we project the above equation in the xy-plane and obtain 
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0, (3.27) 



(3.28) 

From eq. (3.27), we obtain two eigenmodes, to be labeled as the plus "+" mode and minus "— " mode. The refractive indices 
of the two modes are given by 
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2 (r)xx + Vyy) ± V - %y) 2 + 4j ) 
n± = - . (3.29) 

We write the mode polarization vector as E± = E±t + E± z z in the xyz-ir&me, with the transverse part given by 

E *r= (i + j^i/a ^- 1 )- ( 3 - 3 °) 
where 

K ± - ( — \ - nvv - - ~ Vvy ^ - E + 4rix v rt y x (331) 

describes the polarization state of the two eigenmodes. 

From the dielectric tensor of relativistic streaming pair plasma given by eqs. (3.14)-(3.16), we obtain the tensor compo- 
nents in the xyz coordinate system: 

ixx = e x 'x' cos 2 ®b + t z 'z' sin 2 9 B ~ (e x 'z' + sinfls cos9 B , 
e yy ~ e y'y', 

e zz = e x ' x ' sin 2 9 B + e z 'z' cos 2 9 B + (e x > z > + e z > x >) sin 9s cos#s, 

£xj/ = -Cya: = <Vi/' COS 6» S - £ z /y/ sin 9 B , 

£i Z = = (e x ' x i — £ z <z<) sintfs coa9 B + e x ' z ')(cos 2 6>s — sin 2 9 B ), 



r 



V z = —e zy = e y r x > sm9 B + e y > z r cos9 B , (3.32) 



Combining the above equations with eq. (3.27), we find that r\ xy = —r] yx , and r\ yx is almost purely imaginary (except very 
close to cyclotron resonance). We define the polarization parameter, /3 po i, as 

. r) xx — i]y y .e y > y > — e x ' x > cos 2 9 B — t z 'z' sin 2 9 B + (e x > z > + e z > x >) sin 9 B cos 9b „„, 

Ppoi = -1 — ^ -1 —. 5 . Q \ (3.33) 

2?7^ 2(e y / x / cost's - £3,'^' sin 9 S ) 

Then eq. (3.31) can be written as 



iK± = ySpoi T sign(7 7ya; ,,)^ / //32 ol + l. (3.34) 

Here sign(%a;,i) means the sign of the imaginary part of ij yx . Obviously, when \P po \\ 3> 1, the two eigenmodes are linear 
polarized, while for \f3 po \\ = the two modes are circular-polarized. 

Consider a cold pair plasma with AiV = N p — N e <C N and A7 = 7 P — 7 e <§C 7. When the Lorentz-shifted fre- 
quency, fu>(l — f3cos9 B ), is much less than the cyclotron frequency lo c , i.e. for r <JC r cyc or A = w c /[7w(l — f3cos9 B )] = 
u 1/2 7 _1 (1 - (3 cos 6>s) _1 > 1, we have 

(1 - 6»| 7 2 )AA^/iV - A7/7' 

Here we assume 9 B <C 1, so that A ~ 4w 1 / 2 7(l + 9 B r y 2 )~ 1 . After the photon passes through the cyclotron resonance, r 3> r cyc 
or A < 1, the polarization paramerter is given by 

/3p ° 1 ~ 0%1 2 (3 - 0| 7 2 ) A7/7 - (1 - 9%Y)AN/N ■ (3 - 36) 
These expressions are useful for understanding the effect of mode circularization (section 4.3). 



R ~ —™B I \± -r " B 1 ) f o nc\ 

p P o1 - Ti 02 „, 2 \ a at/ at TTZTZr (a.jo) 



3.4 Evolution of Mode Amplitude 

In the xyz frame [with z — k, B — (— sin#s,0, cos 9s) in this frame], we know there are two wave modes: "+" mode and 
"— " mode. It is convenient to introduce a mixing angle, 9 m , via tan# m = l/(iK + ), so that 

tan26» m =/3- 1 1 . (3.37) 
In the xyz frame, the transverse components of the mode eigenvectors are 

E + =( iC ° S <!™), E^f-^Y ( 3.3 8) 



sinf m / v cost 

In the fixed XYZ frame (see Fig. 1), they become 

E _ ( i cos 9 m cos 4> B — sin 9 m sin 4> B \ _ ( — i sin 9 m cos S — cos 9 m sin S 1 , , 

+ \ icos9 m sin <f> B + sin# m cos cj> B /' \ — i sin 9 m sin 4> B + cos 9 m cos (f> B ' 
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The general wave amplitude can be written as 

Z)= a+e++a - e - 

Substitute this into the wave equation, we obtain the mode amplitude evolution equation: 



. d 



A, 
A- 



-Ak/2 + <j)' B sin 20 m i0' m + 4>' B cos 20 m 
-i0' m + 4>' B cos 26 m Ak/2 - <f>' B sin 20 „ 



A- 



(3.40) 



(3.41) 



where the superscript (') specifies d/ds, AA; = k + — k- = Anto/c, and we have subtracted a non-essential unity matrix from 
the above. This equation generalizes the special cases (where only m or 4> B varies) studied in Lai & Ho (2002,2003) and van 
Adelsberg & Lai (2006), and it is useful for understanding the effect of mode coupling (section 4.2). 



4 SOME IMPORTANT PROPAGATION EFFECTS 

With the equations derived in previous sections, we can now identify several key physical effects relevant for the evolution 
of wave polarization. We consider the "weak dispersion" region where the wave frequency is much larger than the plasma 
frequency in the plasma rest frame and the refractive indices of the two natural wave modes are very close to unity. So we do 
not discuss the refraction effect here. The detail discussion about refraction effect can be found in Barnard & Arons (1986). 



4.1 Cyclotron Resonance/ Absorption 

Cyclotron resonance occurs when the wave frequency in the electron/positron rest frame is close to the cyclotron frequency: 

eB 

Co — jlj (1 — /3 cos B ) = lu c = . (4.42) 

mc 

The eigenmodes at cyclotron resonance point are always two circular polarized modes (marked as "+ " for the left-handed 
circular polarized mode and "— " for the right-handed one). Since the electrons and positrons have different directions of 
gyration (one is right-handed, the other one is left-handed) , the right-handed circular polarized mode is absorbed by electrons 
while the left-handed circular polarized mode absorbed by positrons. For right-handed circular polarized mode, the scattering 
cross-section by electrons in the electron rest frame (the physical quantities in the rest frame are marked by " ~ " ) is 

9 e 2 

a- ~ (2tt) 2 8(lj - Wo). (4.43) 

mc 

The opitical depth of this mode in the rest frame is 

f- = J N e ads (4.44) 
Since the optical depth is Lorentz invariant, and 

N e = -y- 1 N e , ds = 7e (l- /3 e cos6) s )ds, (4.45) 
the optical depth in the "lab" frame is 

t- =f- = J N e a{l- /3 e cos0 B )ds. (4.46) 
For a simple model, we set: 

B(r)^B,(^)\ N. { r) (4.47) 
with B„ the surface magnetic field. Thus the optical depth is given by (e.g. Rafikov & Goldreich 2005) 

r- ~ ^(1 - /3 e cosfls)^ ~ 0.62^ 1 1 / 2 3 7 - 1 / 3 z,- 1/3 P 1 ; 1 (l - ft cosfe) 2 / 3 . (4.48) 

with S^ia = B*/(10 12 G), ia, = *V(10 9 Hz). From eq. (4.42) we can find the resonance radius of the electron 
r e , C yc/R* = 1.8 x 10 3 S, 1 / 3 z/- 1/3 7e - 1/3 6'- 2/3 . (4.49) 
The optical depth of the left-handed circular polarized mode caused by the scattering of positrons is similarly given by 

T+ ~^,fr(l-&CO60 B )^. (4.50) 
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with r\ p = N p /Ngj and r PyCyc denned by eq. (4.49) except using 7 P instead of j e . 

When there is an asymmetry between electrons and positrons (different density and/or different 7), the optical depths 
of the two modes are different: 

with t ~ r+ ~ t-. Now consider a linear-polarized photon propagating through the cyclotron resonance region. The mode 
evolution is non-adiabatic (which is always the case since the resonance happens after the polarization limiting radius; see 
sect. 4.2). Before the resonance, the total intensity is 

h = h,+ + Ji,-, with /;,+ = 7i,_ (4.52) 

which means that the intensities of the two circular-polarized modes are the same. The wave intensity after the cyclotron 
absorption is 

If = 7 f ,+ + 7 f ,_ = I u+ e- T + + h,-e- T - . (4.53) 

Because of the difference between r+ and r_, the final intensities of the two circular-polarized modes are different. Thus 
circular polarization can be generated: 

V f = 7 f ,+ - If,_ = Ii,+e- T + - h,-e~ T - . (4.54) 
V f _ e" T + - e- T - 



I f e- T ++e~ T - 



(4.55) 



When T± « 1, V f /I f = -Ar/2 = r _ aat). 

We can also obtain the same result formally by using the Stokes parameters evolution equation (3.25). Since electrons 
and positrons have slightly different 7, the cyclotron absorptions caused by electrons and positrons occur at different radii. 
We analyse them separately. Consider the cyclotron absorption caused by electrons first. Near the resonance, with 

x = — Ve ' cyc , |x|«f, (4.56) 



we have 

ui7 2 (1 -f3cos8 B r 2 ^ ( — — I = (f - xf ~ f - 6x, (4.57) 



' e,cyc 

where we have assumed B oc r -3 . The imaginary part of F\\ and F12 in eq. (3.25) are 

Fl - - - ^(/..") - (-^y = §£r i + ( 3^ d r (4 - 58) 

Fia.i = /e,ia,i =Im(/ e ,i 2 ) ~Im[-/e,n(l + a;) 3 ] ~ -/ e ,u.i- ( 4 -59) 

Also <C |/ e ,n,i| near the resonance, so we neglect it. Thus, the evolution equation for I and V in eq. (3.25) are simplified 
to 

d7/dr = -feo/e,n,i7 + fco/e,n,iV, 

dV/dr = fco/ e ,ii,i/-fco/ e ,ii,iV. (4.60) 
Then we have 

d7+/dr = 0, 

d/_/dr = -2k feM,iI-, (4.61) 

with 7 + = ( J + V) /2 the intensity of left circular polarized mode and 7_ = ( I — V) /2 the right one. The solution of these 
equations is 

If = h+ + h,-e~ T - 

V f = + Ii,-e- T - , (4.62) 

e h, 

( 2fco/ e ,ii,idr = y??e(l-/?eCOs6> S )^^, (4.63) 

»/ acrossCR ' 
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in agreement with eq. (4.48). For the cyclotron absorption by positrons, the analysis is exactly the same, except / Pj i2 = 
+ x) [i ~ / p ,ii. The intensities evolution equations are 

d/+/dr = -2fco/ p ,ii,»i+, 

d/_/dr = 0. (4.64) 

Including both cyclotron absorption by electrons and positrons, the intensity and Stokes V parameters after the resonance 
are 

// = h+e- T + + Ii,-e- T - 

V f = Ii,+e- T + - Ii,-e~ T - (4.65) 
with 

f 1-k r 

r+ = / 2feo/ P ,n, i dr = — Th,(l-A,coBfe)-£^. (4.66) 

Thus, our evolution equations for the mode and Stokes parameters derived in section 3 automatically include the correct 
physics of cyclotron absorption by electrons and positrons. 

4.2 Wave mode coupling 

Wave mode coupling happens near the "polarization limiting radius", r p \, where the mode evolution changes from adiabatic 
to non-adiabatic, i.e., from F a d(r < r p i) > 1 to r at j(r > r p i) < 1. Generally, this is caused by the rotation of the pulsar. 
Obviously the concept of wave mode coupling is relevant for determining the observed polarization only when the wave mode 
is linear polarized, i.e. r p \ < r cir (see section 4.3). In the process of wave mode coupling, the circular polarization will be 
generated. For r < r c i r (so that 9 m — or 7r/2), the mode amplitude evolution equation (3.41) simplifies to 



d f A+ \ ( -Afc/2 i(j>' B \ [ A 



(4.67) 



(4.68) 



l ds y A- J y -i<j>' B Afc/2 J \ A. 
with Ak — Amo/c. The adiabatic parameter is defined as 
Anuj 

1 ad = 

where 

An = sinfls - /n[l - (cos6» s - $sin6> s ) 2 ]}. (4.69) 
When r < r cyc , An ~ 1/2/,, sin 2 6 B ~ — 2v8 B 2 -y~ 3 , so we have 

r ad = 5.6 x lO^Bia^sV^r 1 , (4-70) 
where we have used <j)' B = F^/ric, r\ c = c/fi, and 

„ sin 2 a cos C — sin a cos a sin C cos tb , , „„ . 

= 4. (4.71) 

1 — (cos a cos (" + sin £ sin a cos i/>) 

Obviously, r a( j 3> 1 means adiabatic mode evolution while r at j <§C 1 non-adiabatic. The condition r at j(r = r p i) = 1 then gives 
the polarization limiting radius 

r pl /R = 1.8 x lOy/^V^B^V'l^r 173 - (4.72) 
Compare r p i with r cy c, we have 

r P ,/r cyc = r, 1/3 7 - 2/3 |F0r 1/3 = 0.215^ /3 ^- 2/3 1 I " 1/3 £ 1- (4-73) 

So in the typcal parameter region (n = 100, 7 = 100, = a few), wave mode coupling always occurs before cyclotron 
absorption. 

To understand the wave mode coupling around r p i, we write 

r ad = x~ n , (4.74) 

with x = r/r p \. According to eq. (4.70), the power-law index n ~ 3 (not exactly 3 because 9 B also varies as r changes). Then 
eq. (4.67) can be simplified to 

IM ( Z", , (4.T5) 



di I i_ I \ -sign(0' s )i -x n J \ A 
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Figure 2. Evolution of the radiation mode amplitudes (top panel) and Stokes parameters (bottom three panels) with power-law index 
n = 3 [see eq. (4.74)]. The solid lines are for A = 1.0 and the dashed lines for A = 0.1. The polarization limiting radius is at x = 1. The 
initial values (at a small x = x\) are = 1, A_ =0, Q = I = I, U = and V = 0. When x <J 0.5, the modes evolve adiabatically. At 
r ~ r p i (or x = 1) the modes begin to couple, generating circular polarization. At d > 1, the Stokes parameters are "frozen". 



where 

A = r pl 4>' B = O.aSsignl^)^ 1 / 3 ^ 1 ^ 3 ^ 1 /^- 2 /^- 1 ^- 1 !^! 273 - (4.76) 

Similar equation is given by van Adelsberg & Lai (2006), except that in their paper the dispersion relation of X-ray is 
dominated by QED effect so that An > 0, while in our case plasma effect dominates the radio wave propagation with An < 0. 
Figure 2 shows two examples of mode evolution with A = 0.1 and A = 1.0, both for n = 3. The photon is 100% linear 
polarization before the wave mode coupling (Here we set it to be O-mode initially). After the wave mode coupling (x 3> 1), 
the polarization states are frozen. In this process, circular polarization is produced. It is obvious that the larger A is, the more 
circular polarization will be generated. Figure 3 shows how the value of A affects the final circular polarization when 
the power index n = 3. For n = 3 and A < 0.1, the final circular polarization is given by the expression 

y// = 2.2sign(^)|A| 3/2 . (4.77) 

For A <; 1, the circular polarization l^j/i is close to 1. Equation (4.77) also shows the relationship between the sign of the 
circular polarization and 4>'b- An increasing <j>B (or (j>' B > 0) corresponds to positive circular polarization while decreasing cf>B 
to negative one. 



4.3 Circularization 

Circularization happens when |/3 po i| ~ 1, and we can define the radius of circularization r c i r by |/3 po i(r = r c j r )| = 1. For 
t 3> 7" c i r , the normal modes become circular-polarized. 

According to eq. (3.35), if r c i r <C r cyc (or A> 1, before cyclotron resonance) and 6bJ 1, the polarization parameter 

1^1 " ^AivA-A 7 / 7 (4 ' 78) 
which means the two wave modes are always linear polarized. However, if r C ir <^ r c yc and #^7 <§C 1, 

\ n2 2 

I' 3 - 1 " ^ A7V/7V A7/7 ' ^ 
so circularization could happen when 




Figure 3. The final circular polarization fraction |V|// after wave mode coupling as a function of A with the power-law index n = 3. 
The linear polarization fraction before wave mode coupling is assumed to be 100%. The dash line depicts the fitting formula [eq. (4.77)] 
for A < 0.1. 



6»b7 = \J (AN/N - A 7 / 7 )/A. (4.80) 

Which implies a very small 9b- This condition could be satisfied when the photon ray is nearly aligned with the magenetic 
field, or when the photon is generated inside the I/7 cone of the radiation beam. 

If the circularization happens after the cyclotron resonance (r c i r 3> r cyc or A <C 1), according to eq. (3.36), the radius of 
circularization is given by 

r cir /R, = 2.2 x 10 3 Bl / 1 lu- 1/3 6- 4/3 -y- 1/3 {AN/N - A 7 / 7 )~ 1/3 . (4.81) 

Here we used B ~ B,(r/i?*) -3 and assumed 9bJ S> 1- The ratio of r c i r and the cyclotron resonance radius [see eq. (4.49)] is 

r cil /r cyc = 1.2 9~ 2/3 (AN/N - A 7 / 7 )~ 1/3 > 1. (4.82) 

Obviously, in the parameter regions we interested in, the circularization radius is typically larger than the cyclotron resonance 
radius and the polarization limiting radius. Thus, this effect does not change the photon polarization state at all. 



4.4 Quasi- Tangential Propagation Effect 

In their study of the X-ray polarization signals from magnetized neutron stars, Wang & Lai (2009) found that as the X-ray 
photon travels through the magnetosphere, it may cross the region where its wave vector is aligned or nearly aligned with 
the magnetic field (i.e., 8b is zero or small). In such a Quasi-Tangential region (QT region), the azimuthal angle of magnetic 
field <f>B changes quickly, the two photon modes (|j and T modes) become (nearly) identical, and mode coupling may occur, 
thereby affecting the polarization alignment. This Quasi-tangential Effect generally happens at a few J?„ for surface X-ray 
emission. The physical mechanism is similar to the wave mode coupling effect discussed in section 4.2 [see the mode evolution 
equation (2.11) in Wang & Lai (2009)], except that the magnetic field plays an important role. 

In the radio case, we assume the photon is emitted in the tangential direction of the magnetic field line at the emission 
point (~ 507?,). If the NS is non-rotating, then the k — B angle 9b (= at the emission point) will increase monotonically and 
no QT effect will occur. However, when we consider the rotation of the NS, for some special photons (for example, those with 
small impact angle \ and special ^i) 9b could attain its minimum value at a large radius. As an example, the two bottom 
panels of Figure 6 and 7 (to be discussed in detail in section 5) shows the evolution of 9b, 4>b along the ray for \ ~ 0.5°. We 
see that 9b reaches its minimum value at about s — 7007?, away from the emission point. The azimuthal angle <f>B changes 
very quickly at this radius. The two linear modes strongly couple with each other. The final polarization state after crossing 
this QT region is complicated: </!>pa can be modified significantly and different sign of circular polarization can be generated 
for different geometry, which is different from wave mode coupling effect discussed in section 4.2. In general, the QT effect 
strongly influence the polarization phase profiles when impact angle is very small (see section 5.3). In our case, the QT effect 
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is always coupled with wave mode coupling effect (occurring at almost same place), and the numerical ray integration is 
necessary to account for these effects accurately (see section 5). 

5 NUMERICAL RESULTS 

In section 4 we have discussed various key physical effects related to wave propagation through the magnetosphere. However, 
in many cases these different effects are coupled and not easy to separate. Thus, to produced the observed polarization profiles, 
it is necessary to use the numerical ray integrations to calculate the final wave polarization states. 

5.1 Single Ray evolution 

It is generally accepted that pulsar radio emission is emitted from the open field line region at a few to tens of NS radii (e.g. 
Cordes 1978; Blaskiewicz et al. 1991; Kramer et al. 1997; Kijak & Gil 2003). In this paper, we choose the emission height 
f'cm = 50i?* and assume that at the emission point, the photon is polarized in the k-B plane (or the O-mode, as in the case 
of curvature radiation) , and propagates along the tangential direction of the local magnetic field line (here we do not consider 
the emission cone of angle I/7). For a given emission height r cm , the pulsar rotation phase the direction of line of sight 
£ (which is the k-fl angle), the surface magnetic field B„, and the plasma properties (plasma density parameter rj, Lorentz 
factor of the streaming plasma 7), we can calculate the dielectric tensor at each point along the photon ray, and integrate the 
wave evolution equation (2.10) from the emission point to a large radius (generally we choose ri c /2), beyond the polarization 
limiting radius r p i and cyclotron resonance radius r cyc , to determine the final polarization state of the photon. 

5.1.1 Symmetric pair plasma 

We first consider the case of symmetric pair plasmas, i.e., the electrons and positrons have the same Lorentz factors (7^ = 7 e , 
or A7/7 = 0) and densities (AN/N — 0). In this case, the eigenmodes are always linear polarized (mixing angle 9 m = 0° or 
90°). Figure 4 shows an example of the photon polarization evolution along its trajectory. We can clearly find the wave mode 
coupling effect (at r p \ ~ 800 R*) and cyclotron absorption effect (at r cyc more than 1000i?*). The final polarization position 
angle 0pa is determined by </>s(r p i) [see eq. (5.87)]. It is obvious that near the polarization limiting radius, r a( j oc (r/r p i)~ n 
and n ~ 3, so that as discussed in section 4.2, the final circular polarization is determined by the value of A [see eq. (4.77)]. 
Since we are dealing with a symmetric pair plasma here, cyclotron absorptions do not change the polarization state (but 
decrease the total intensity) . 

Figure 5 give some other examples of the evolution of Stokes parameters with different plasma density rj and Lorentz 
factor 7. Different rj and 7 correspond to different r p i [according to eq. (4.72), lower 77 and higher 7 corresponds to a smaller 
r p i] , so that the final 0pa is different too. In all the above cases, the final polarization state changes significantly compared to 
the original state, not only the linear position angle but also the circular polarization. 

At the special parameter region of the initial rotation phase Vl/i, the QT effect (see section 4.4) can strongly affect the 
final polarization state. Figure 6 shows the photon evolution for = —9° (the other parameters are the sames as in Fig. 4, 
e.g. the impact angle \ — 2°). Note that in contrast to Fig. 4, here the k — B angle 6 b does not vary monotonically along 
the ray. There exists a QT region around s ~ 700i?, , where b is minimum and <f>B is changing very quickly. As discussed in 
section 4.4, the final 4>p\ and circular polarization are different from the prediction of pure wave mode coupling effect (which 
is the case in Fig. 4 where QT effect does not occur). For the photon evolution with a smaller photon impact angle \ = 0-5° 
(but the initial rotation phase and other parameters are the same as in Fig. 6), the QT effect is stronger, as shown in Figure 7. 
Note that even the sign of the final circular polarization in this figure is positive, as a result of the strong QT effect. 

5.1.2 Asymmetric pair plasma 

If the electrons and positrons of the magnetospheric plasma have different velocities and/or densities, the wave eigenmodes 
cannot always be linearly polarized. As discussed in section 4.3, before cyclotron resonance the natural modes are linearly 
polarized [see eq. (4.78)] for 6bJ >• 1. After the cyclotron resonance, the natural modes become elliptical polarized. In section 
4.3, we have defined a circularization radius r c i r where the polarization parameter |/3 po i| = 1 [see eqs. (4.81) and (4.82)]. For 
r ^> r c i r , the natural modes becomes circular polarized. 

According to eq. (4.73), for typical plasma parameters of interest in this paper, rj ^> 100 and 7 i; 100, wave mode coupling 
always occurs before the cyclotron resonance (r p i < r cy c)- Thus, circularization always happens after wave mode coupling, at 
which point the wave polarization state is already frozen. Therefore the change of natural mode does not affect the observed 
polarization state. Figure 8 shows the photon evolution in an asymmetric pair plasma. Note that the mode mixing angle 6 m 
changes from 0°/90° to 45° after the cyclotron resonance, but the polarization state does not change since r c i r > r p \. 
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s/R* s/R, 

Figure 4. An typical single photon evolution across the magnetosphere. The horizontal axis s/ R* is the photon distance away from the 
emission point. On the left panels, m is the mode mixing angle defined by eq. (3.37), r ad is the adiabatic parameter defined by eq. (4.70) 
(the dashed line is the power-law fit of r ad around r p i, which is r ad <x s -2 39 ), Ax and Ay are the wave amplitudes in the fixed XY Z 
frame. A+ and A_ are the mode amplitudes, and 9g is the angle between k and B. On the right panels, 7, Q, U, V are the Stokes 
parameters, ipPA = 0.5tan~ 1 (C//Q) is the linear polarization position angle (solid line) and 4>b is the azimuthal angle of the B field 
(dashed line). The initial polarization is assumed to be in the ordinary mode, with A+ = 1, A_ = 0. For this example, the parameters 
are: surface magnetic field B, = 10 12 G, NS spin period P = 1 s, wave frequency v = 1 GHz, plasma density parameter r\ = N/Nqj = 400 
(N = N e + N p and N e = N p ), Lorentz factor 7 = 100 (with A7/7 ~ 0), inclination angle a = 30°, impact angle \ = 2°, initial rotation 
phase \Pi = 0°, and emission height rem = 50i?*. 



von Hoensbroech et al. (1998) studied wave modes in a pure electron plasma. They assumed that the background plasma 
has a much lower Lorentz factor (e.g. 7b g = 1.7) than the Lorentz factor of the radiating beam. In this case r c ir may be close 
to r p i and circular polarization may be generated around r c i r . Note that they did not calculate r p i but simply assumed that 
the final photon polarization is determined by the normal mode at some fixed r p i. 



5.2 Polarization Profiles of Pulsar Emission Beam 

Having understood the main features of polarization evolution along a single ray, we now proceed to calculate the polarization 
profiles of pulsar emission beam. To do this, one needs to know the emission height as a function of pulsar rotation phase. 
For simplicity, in this paper, we assume that all emissions are from the same height, at r cm = 507?», and defer the results 
for emissions from a range of heights to a future paper. For a given emission height r cm , the pulsar rotation phase "Ji, the 
inclination angle a and the direction of line of sight £ (which is the k-tt angle), we can find the position of the emission point 
ri where the tangential magnetic field line direction is along the line of sight. This emission point, ri = (r cm , (9,i, (j> r i), is 
given by (in the fixed XY Z frame) 

6 r \ = -y - ^sin" 1 (j sin6 Vi) , 4>n = 4>^i, (5.83) 

where (9^, (fr^i) is the initial direction of the dipole magnetic momentum \i i and can be found in eq. (2.5) (with ^ given by 
^i). We consider emissions only from the open field line region, i.e., the angle between r\ and /x ; should be less than sj r eul /ri c . 

For given r om , a, f, and initial polarization state (ordinary mode), we determine ri and calculate the final observed 
Stokes parameters by integrating along the ray. When the phase i&j varies due to NS rotation, we can observe photons 
from different emission points and the final observed Stokes parameters will change with the rotation phase — this is the 
pulsar polarization profile. If we neglect the propagation effect, the observed position angle <^pa can be described by the 
Rotating- Vector-Mode (RVM) (see Radhakrishnan & Cooke 1969) as 

, ,-r x -1 — sin a sin 'I'; , 

0pa = <?V = M^i) = tan r^ 2 ST- 5 - 84 

sin Q cos a — cos Q sin a cos W; 
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Figure 5. Same as the right panels of Fig. 4, except for different plasma density r\ and Lorentz factor of the streaming motion 7: the 
solid lines are for r\ = 400, 7 = 100 (same as Fig. 4), the dashed lines for 77 = 400, 7 = 300 and the dotted lines for r\ = 100, 7 = 100. The 
thin line in the bottom panel is for <j>B- It is obvious that for lower density and/or higher 7, the wave mode coupling occurs at smaller 
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Figure 6. Same as Fig. 4, except for a different initial rotation phase \Pj = —9°. The impact angle is the same as in Fig. 4, \ = 2°. 
Note that in contrast to Fig. 4, here the k— B angle 9b does not vary monotonically along the ray. There exists a quasi-tangential (QT) 
region around s ~ 700i?», where 9b is minimum and 4>b changes very quickly. Strong circular polarization is generated here and the 
final PA angle <j>VA can not be predicted by simply using eq. (5.84). 
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Figure 7. Same as Fig. 6, except for the impact angle \ = 0.5°. The very smaller impact angle makes 4>b change more quickly around 
s ~ 700-R,, so that the QT effect is stronger than the case shown in Fig. 6. Note that the sign of final circular polarization is different 
from the cases in Fig. 4 and 6. 



The basic assumption of the RVM is that the radiation is emitted with polarization in the plane of the field line curvature 
(i.e. the k-B plane) and this polarization direction is unchanged during the propagation. However, as seen in section 5.1, the 
final polarization state can be modified compared to the initial one because of the propagation effect in the magnetosphere, 
so that the final PA profile can deviate significantly from the RVM model. 

Figure 9 shows a typical example of the phase evolution of the intensity and polarization, taking into account of all the 
propagation effects. The total intensity is only affected by cyclotron absorption, and a higher plasma density leads to stronger 
absorption. We see that the relative intensity I /Io varies with the rotation phase \&i , simply because the wave passes through 
different paths in the magnetosphere for different ^i. For illustrative purpose, we consider the initial intensity profile Io, given 
by a Gaussian centered at \&i = 0: 



/o(*i) = exp(-4^Irr2*?/(*Lx)). 

Here *& ma x is the initial phase of the photon from the edge of the open field region and is given by 
cos # opcn — cos £ cos a 



cos $ n 



sin £ sin a 



(5.85) 



where opan ~ r cm /n c is the half-cone angle of the open field region at emission height r om (here we simply assume the open 
field region is always the same as the n || SI case). Since I /Io depends asymmetrically on ^i, the observed intensity / is no 
longer a Gaussion. Non-gaussion profiles have been observed in many pulsars, and the phase-dependent cyclotron absorption 
illsutrated here is a possible explanation. 

The final polarization profiles are also strongly affected by the propagation effects. When the plasma density is not so 
high, and/or the impact angle \ ls n °t so small [compared to the half cone angle of the emission beam from the open field 



region 6» bcam ; e.g., in Fig. 9(a), \ = 5° while bc 



l-59 opcn = 8.8°], the wave mode coupling effect is not strong and the final 



circular polarization is not very high. In this case, the final linear polarization position angle is determined by the azimuthal 
angle of B field at the polarization limiting radius <j>B(r p i): 



B (r p i) ~ 7t + (j>n(r p i) = 7t + <^ M (*i + r p \/r\ c ). 



(5.87) 



Here we have used the approximation of B(r) = — /x/r 3 since r p i 3> R*- In general, the polarization limiting radius, r p i, does 
not vary significantly with different rotation phase ^i. Thus the final PA profile just shifts by the amount 



rpi/ric ^ 0.1 



Go 3 ) 



UO 12 GJ UGHzJ 



1/3 



(5.88) 
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Figure 8. Same as Fig. 4, except for an asymmetric pair plasma with A7/7 = 0.2, AN/N = 0. Since the Lorentz factors of electrons 
and positrons are different, their cyclotron resonance occur at different radii. Note that after the cyclotron resonance, the natural modes 
become circular polarized (9 m = 45°). 



compared to the RVM model. The final circular polarization is always single sign in this case. We can easily find the relationship 
between 0pa and the sign of circular polarization. According to eq. (5.87), the monotonicity of the final PA angle ^>pA(^i) is 
given by 

_ d^ s (r pl ) d* _ d(j> B {r pl ) 

d*i d* d*i ~~ d* ' 1 ' ' 

Here 9 = + s/r\ c [see eq. (2.2)], so that d^/d^i = 1. The sign of the final circular polarization (generated by wave mode 
coupling effect) is determined by d<^s(r p i)/ds [<j>' B in eq. (4.77)] and: 

d<M r pi) _ d^ B (r p i) d£ = d<ft B (r p i) J_ 
ds ds d* He' 

So that d(/>pA/d*l/i always has the same sign as d(j>B(r p \)/ds. According to eqs. (4.77) and (5.90), we can find that because 
of the wave mode coupling effect, a monotonically increasing <f>pA leads to a positive V while a monotonically decreasing <j>PA 
gives a negative V. This relationship has been observed in some conal-douple type pulsars; see section 5.3. 

The polarization profiles can also be quite different from the RVM prediction, especially in the case of low impact angle 
and/or high plasma density. Figure 9(b) give some examples for the impact angle \ = 2°. For the low density case of r\ = 10, 
the final PA profile can still be approximated by a simple shift from the RVM model. However, for higher density [r\ — 200, 
400), the final PA profile is not just a simple shift compared with the RVM model. For example, the PA profile of r\ = 400 
case has a 90° jump within 1° around ^ ~ —9°, where the linear polarization L/I is close to while the circular polarization 
|V|/7 reaches almost 100%. In this region the QT effect (discussed in section 4.4) plays an important role in determining the 
final polarization state (see Fig. 6). 



5.3 Two-Dimensional Polarization Maps of Pulsar Emission Beam 

For a given pulsar, observation at different line of sight (i.e., different ( or impact angle \) would obviously result in different 
intensity and polarization profiles. Figure 10 gives an example of the two-dimensional polarization map of the observed Stokes 
parameters, produced by varying \ and while keeping all other parameters fixed. As discussed before, the final total 
intensity I is only affected by cyclotron absorption, while the linear and circular polarizations are modified by wave mode 
coupling effect and QT effect, and can deviate significantly from the prediction of RVM model. Figure 11 shows four profiles 
with four different impact angle \i corresponding to four sections of Fig. 10. These four sections represent three typical final 
polarization states produced by the propagation effects: 

(i) For a relatively large impact angle |x|, the final PA profile can be obtained by a small shift from the RVM profile, 
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Figure 9. The intensity and polarization profiles computed by ray integrations. We use two different impact angles: (a) x = 5° on the 
left panels, and (b) \ = 2° on the right panels. The solid lines arc for plasma density parameter r) = 10, the dashed lines for r) = 100 
and the dot-dashed line for r\ = 400. The top panels show the total intensity profiles, where we have adopted (for illustrative purpose) a 
Gaussian initial intensity profile /rj(*i) = exp(— 4\/ln 2 l I' 2 /*I' max )) (this initial profile is shown shown as dotted lines, almost coincident 
with the solid lines). The second panels form the top show the modification factor I /Iq due to propagation effects. The bottom three 
panels show the linear polarization fraction L/I, circular polarization V/I and the position angle of the LP, 0pa- In the bottom panels, 
the dotted lines (almost coincident with the solid line on the left panel) show the prediction from the RVM model: 4>PA = 4>ni — 4>Bi- In 
these calculations, the initial polarization states are all ordinary mode, and the other parameters are: surface magnetic field B t = 10 12 G, 
NS spin period P = Is, wave frequency v = 1 GHz, Lorentz factor 7 = 100 (with A7/7 ~ 0), inclination angle a = 30°, and emission 
height r cm = 50-R* . 



of the amount r p i/r\ c [see eqs. (5.87) and (5.88)]. Figure 9 and the x = — 5° column of Fig. 11 depict some examples. The 
final circular polarization is always of single sign: a monotonically increasing 4>pa leads to a positive V while a monotonically 
decreasing </>pa gives a negative V [see eqs. (4.77) and (5.90)]. This behavior is consistent with observations of the double 
cone emission of some pulsars ( "conal-double type pulsars"), where a correlation between the sense of CP and the sense of 
PA variation was found (see Han et al. 1998). 

(ii) For relatively small impact angle, the final PA profile is very different from the RVM prediction - the middle two 
columns (x = -1.9° and x = -1°) of Fig. 11) give some examples. It is clear that there always exists a special line of sight 
(x — Xjump), for which the PA profile has a 90° jump (where Q — while U changes signs [see the \ — —1.9° column of 
Fig. 11]. The large jumps in V/I and <^>pa are caused by the QT effect. For |x| < Xjump, the PA is not necessarily an monotonic 
function of *5>i. Nevertheless, the final CP retains a single sign, which is the same as the case with large jx|- 

(iii) For very small impact angle (|x| <C Xjump), the QT effect is much stronger, so that the final PA profile is very different 
from the prediction of RVM model and even the CP does not stay at a single sign (see the right-most column of Fig. 11). 

The above three types of polarization behaviours always exist for different pulsar and plasma parameters (e.g., B, , Po, 
v, r\, 7, a and r om ). Different parameters just modify the position of Xjump and the initial rotation phase where the 90° jump 
in PA occurs, while the basic morphology of the emission beam does not change. 



6 CONCLUSION AND DISCUSSION 

We have studied the evolution of radio emission polarization in a rotating pulsar magnetosphere filled with relativistic stream- 
ing pair plasma. We quantify and compare the relative importance of several key propagation effects that can influence the 
observed radio polarization signals, including wave mode coupling, cyclotron absorption, propagation through quasi-tangential 
(QT) region, and mode circularization (due to asymmetric distributions of electrons and positrons). We use numerical inte- 
gration of the photon polarization along the ray to incorporate all these propagation effects self-consistently within a single 
framework. We find that, for typical parameters of the magnetospheric plasma produced by pair cascade, and for an initially 
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Figure 10. Two-dimensional polarization map of pulsar emission beam. The four panels correspond to the four Stokes parameters /, 
Q/I, U/I, V/I. The filled region is the emission beam from the open field region. The values of the four Stokes parameters are shown as 
different colors (with red for positive value while blue for negative). The map is obtained by computing the observed wave polarization 
for different impact angle \ (which varies from — Sbeam to Sbeamj here <9beam — 8.8°) and the rotation phase The other (fixed) 
parameters are: surface mag netic field B* = 10 12 G, NS period P = 1 s, wave frequency v = 1 GHz, plasma density r) = N/Nqj = 400 
(N = N e + N p and N e = N p ), Lorentz factor 7 = 100 (with A7/7 ~ 0), inclination angle a = 30°, and emission height r cm = 50-R*. 
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Figure 11. Intensity and polarization profiles of pulsar emission beam. The four columns correspond to four fixed impact angles \ = — 5°, 
— 1.9°, —1° and —0.5°, respectively, which are four sections in Fig. 10. In each column, we plot the Stokes parameters I/Iq (top panel), 
L/I (second panel), V/I (third panel), </>pa and (j>Bi (bottom panel, the solid line for 4>PA and the dashed for ipBi)- Note that there exists 
a 90° jump of <j)pA near = —9° in the bottom panel of the \ = —1.9° column. The other parameters are the same as in Fig. 10. 



100% linear polarized radio wave, the final intensity and polarization position angle are modified by the propagation effects, 
and significant circular polarization can be generated. 

We find that the most important propagation effects are cyclotron absorption, wave mode coupling and QT effect. 
Generally, cyclotron absorption occurs after the wave mode coupling [r cyc > r p i, see eq. (4.73)]. Thus, it only changes the 
total wave intensity and does not modify the wave polarization (4>pa, V/I). For a large impact angle |x| and/or relatively 
low plasma density, the final wave polarization angle (j>pA is determined by the azimuthal angle of B field at the polarization 
limiting radius r p \, and the observed circular polarization is determined by the value of A = r p i(f)' B at r = r p i [see eqs. (4.76) 
- (4.77)]. In this case, the observed 4>pa profile is similar to the prediction of the Rotating- Vector Model (RVM), except for 
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Figure 12. Intensity and polarization profiles of PSR J1920+2650. In the lower panel, the upper thick line is the total intensity (/), 
the dash-dot-dashed line is the linear polarized intensity L, and the thin line is the circularly polarized intensity V. This pulsar shows 
extremely strong circular polarization with V/I ~ 64% in its first major component. 



a phase shift by the amount r p \/r\ c [see eq. (5.88)]; the circular polarization has a single sign across the emission beam (see 
Fig. 9). For a small impact angle and/or high plasma density, the QT effect becomes important, the final polarization profiles 
are more irregular: a 90° sudden jump in PA may occur at certain phase, accompanied by large circular polarization. For very 
small the circular polarization may change signs for at different phases (see the right column of Fig. 11). 

In this paper, we have adopted the simplest (and minimum) assumptions about the property of the magnetospheric plasma 
and the intrinsic radio emission mechanism (see below). Nevertheless, our results already show great promise in explaining a 
number of otherwise puzzling observations: 

(i) It has been observed that in some single-pulse pulsars, the intensity profile deviates from the single guassion shape. One 
possible reason is that cyclotron absorption depends on the rotation phase (because the ray passes through different region 
of the magnetosphere) , as discussed in section 5.2. Thus, even when the initial intensity profile from the emission beam is a 
Gaussion, the observed profile can be non-gaussion. 

(ii) For the so-called conal double type pulsars, which in our model corresponds to large impact angle \, the relationship 
between the single sign of the circular polarization and the derivation of </!>pa (see Han et al. 1998) can be easily understood 
by the wave mode coupling effect. According to eqs. (4.77) and (5.90), an increasing 0pa corresponds to the left-hand circular 
polarization (V > 0) while a decreasing cSpa corresponds to the right-hand (V < 0) one . 

(iii) According to our calculation, there exists a special impact angle Xjump, where the observed <^>pa profile has a 90° jump 
(orthogonal polarization mode) and this is accompanied by the maximum circular polarization. This feature may be helpful 
to explain the polarization profile of PSR J1920+2650 (Figure 12; see Han et al. 2009). (iv) For a very small impact angle, 
which corresponds to the core emission, the QT effect can cause the sign reversal of circular polarization, which is observed 
in the core components of many pulsars (e.g. Radhakrishnan & Rankin 1990; Han et al. 1998; You et al. 2006). 

Our calculations in this paper have relied on several simplifying assumptions. For example, we have assumed that the radio 
emission is from the same height for different rotation phases, that the density parameter (r) = N/Nqj) of the magnetospheric 
plasma is constant everywhere in the emission cone and along the photon trajectory, and that the plasma electrons and 
positrons have the same for bulk Lorentz factors. In future works, we plan to consider models with varying emission heights, 
as well as non-trivial electron/positron spatial and velocity distributions. We did not include the small but finite emission 
cone (angle I/7) in our model, and assumed that the initial polarization of photon is always O-mode for different rotation 
phases. However, different emission mechanisms could give different initial polarization states. We will also be interested in 
studying the propagation effects on the individual/subpluse emissions, since they may more directly reflect the underlying 
radio emission processes. 



1 Here we define the circular polarization as seem from the receiver, which is different from the defination in electrical engineering used 
by Han et al. (1998) 
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